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ABSTRACT 


This report describes a series of low speed airfoil designs based 
on modifications to the NACA 64^-212 airfoil. Designs are based on 
potential flow theory. This report describes one of a series of airfoil 
modifications carried out under Contract NAS 2-8599, Application of 
Multivariable Search Techniques to Optimal Wing Design in Non-Linear 
Flow Fields. Mr. Raymond Hicks of National Aeronautics and Space 
Administration's Aeronautical Division, Ames Research Center, served 
As contract monitor for the present study. 


AN INVESTIGATION ON THE EFFECT OF 


SECOND-ORDER ADDITIONAL THICKNESS DISTRIBUTIONS TO 
THE UPPER SURFACE OF AN NACA 6^-212 AIRFOIL 

by Donald S. Hague and Antony W. Merz 
Aerophysics Research Corporation 

SUMMARY 

An investigation has been conducted on the Lawrence Radiation Center, 
Berkeley, CDC 7600 digital computer to determine the effects of additional 
thickness distributions to the upper surface of an NACA 64^-212 airfoil. 
Additional thickness distributions employed were in the form of two second- 
order polynomial arcs which have a specified thickness, y, at a given 
chordwise location, x. The forward arc disappears at the airfoil leading 
edge, the aft arc disappears at the airfoil trailing edge. At the 
juncture of the two arcs, x = x, continuity of slope is maintained. 

The effect of varying the maximum additional thickness and its chordwise 
location on airfoil lift coeeficient, pitching moment, and pressure dis- 
tribution was investigated. ; Results were obtained at a Mach number of 
0.2 with an angle-of-attack of 6° on the basic NACA 64^-212 airfoil. All 
calculations employ the full potential flow equations for two dimensional 
flow. The relaxation method of Jameson is employed for solution of the 
potential flow equations. 

Increases in the rearward location of the maximum additional thickness 
and increases in the magnitude of the additional thickness both produce 
increases in the airfoil lift coefficient. Conversely moving the location 
of maximum thickness forward or decreasing the maximum thickness both 
reduce the magnitude of the quarter chord pitching moment. The magnitude 
of the largest pressure peak varies in a complicated manner with maximum 
additional thickness and its chordwise location. For maximum thickness 
locations forward of the 2/3 chord additional thickness initially produces 

•7 

a reduction in pressure peak with a lift coefficient increase. With 
larger amounts of additional thickness pressure peak value and lift 


coefficient rise together. For maximum thickness locations aft of the 
2/3 chord location additional thickness produces a monotonic rise in both 
lift coefficient and pressure peak magnitude. A consequence of this 
behavior is that for a given lift coefficient value the peak pressure can 
be minimized by careful selection of the location of maximum thickness 
and its magnitude. Generally as the lift coefficient rises the maximum 
thickness location moves aft. For a of 1.2 the optimal location for 

maximum thickness is at the quarter chord. For a C. of 1.8 the optimal 

fi L 

.location is approximately at the half chord. 

It should be noted that viscous effects are neglected in the present 
analysis. At the higher lift coefficients the effect of viscosity could 
be significant. Further investigations incorporating a viscous flow 
model are therefore desirable. 

INTRODUCTION 

f; 

The National Aeronautics And Space Administration and others are 
currently conducting a series of theoretical and experimental studies 
to define airfoil sections having improved performance from the aspects 
of lift, drag, pitching moment or pressure distribution characteristics, 
refs 1 and 2. Analytic investigations ysing airfoil surface repre- 
sentations based ou high-order polynomials may result in impractical 
profiles, for example, very thin trailing edge thickness distributions 
or severe reflexes in the profile. The present study employs low-order 
polynomial arcs of second-order whose characteristics are selected to avoid 
such problems. Optimization studies using multivariable search techniques, 
reference 3, generally indicate that shape changes which provide increased 
lift produce unfavorable changes in moment characteristics. Conversely 
profile changes which improve the moment characteristics decrease the lift 
coefficient. With the low-order model of the present investigation a 
systematic examination on the effect of .profile changes can be carried 
out. This was accomplished and the trends revealed by optimization 
studies were confirmed. An interesting by product of the systematic 



il 

investigation of profile changes is that a gain in lift coefficient can be 
produced while reducing the peak negative pressures. This tends to 
decrease the pressure gradient and hence holds promise for the development 
of practical single component high lift coefficient airfoils. 


MATHEMATICAL MODELS 

Potential Flow Equation 

Potential flow analysis is based on solution of the two-dimensional 
potential flow equation 

(a 2 -u 2 ) 0 xx + (a 2 -v 2 ) 0 yy -2uy 0 xy =0 

where 0 is the velocity potential, u and v are the velocity components 

u - 0 X , t y 

and a is the local speed of sound determined from the energy equation and 
the stagnation speed of sound 

. a o 2 - (tl) (u 2 - v 2 ) 

Solutions are obtained by Jameson's finite difference scheme, reference 4. 


AIRFOIL PROFILE REPRESENTATION 


Basic Airfoil 

Ordinates for the basic NACA 64 1 -212 airfoil were approximated by 

four cubic chain polynomials in the manner of Hicks 

y. = a F 1 + a. x + a 9 x 2 + a, x 3 ; j = 1,2, 3,4 
: 0.1 1. 2. 3 j 

Cat'ii* 'icients in the four polynomial arcs are selected on the following 









X ~ 1 

i = 2 


1 = 3 
1 - 4 


Ax*c represents 
Arc represents 


forward portion of upper surface 

h - V 7 

aft portion of upper surface 


Arc represents 
Arc represents 


forward portion of lower surface 

F 3 = V* 

aft portion 


p 4 = 1 

The coefficients a, are determined by introducing four boundary conditions 

j 

on the airfoil profile in each of the four airfoil arcs. Crout's method for 
triangularisation and back substitution of the resulting systems of linear 


simultaneous equations. Note that if four points are specified on the 
aft portion (i = 2 ok 4), a discontinuity in slope occurs where the poly- 
nomials join. This produces a small ripple in the pressure distribution 
at the juncture point. However, since the juncture occurs at a region 
of small slope (x = .5) the effect is not significant. The approximate 

NAGA 64-212 airfoil developed by this method is presented in Figure 1. 

y 

Additional Thickness 

In the present study additional thickness is limited to the upper 
airfoil surface. The additional thickness has the form 




Ay 00 = 9 

[■•(f)] 

r , - v 2-1 

Ay 00 = y 

[‘ -(H) ] 


These functions are of second-order ranging parabolically with C = |x-xj. 
Additional thickness is zero at the ending edge (x=0) and trailing edge (x=l) 
and has a maximum of Ay=y at x=x. Additional thickness and slope ofs the 
additional thickness are continuous throughout .the interval o<n<l. Second 
derivative of the additional thickness distribution is constant in the 
forward and aft airfoil arcs but has a discontinuity at the -arc junction, 
x=x. It follows that a continuous polynomial representation of the additional 
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thickness distributions, valid in the interval o<n<l would be in the 
form of an infinite series. This type of additional thickness distribution 
is referred to as a "biquadratic" function in recognition of the above 
characteristics. A sequence of biquadratic arcs having varying maximum 
thickness positions are presented in Figure 2. 


PRESENTATION OF RESULTS 
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Additional Thickness Optimization 


Lift Coefficient Maximization 

Maximization of lift coefficient has the form 

0 = Max[C L ] 


!? 


S) 

i 


where 

Si' 
j) 

and the integration is around the airfoil contour. Since, the airfoil 
contour is completely described in terms of the two parameters x and y 

0 = Max [C L ] * Max [C L (x, y)] 

where 

V *M - 


®Ap(x)dx 


m 

tf 

This two variable multivariable search problem was solved by a combination 
of directed random-ray and pattern searches. Ref. 3. Table I presents the 
results of 30 iterations using these search procedures. Lift gains are 
produced at 27 of the 30 iterations and continue to be made at the compu- 
tation termination. 

Optimization has moved the position of maximum thickness to the most 
rearward position allowed, x = 0.9. At termination, lift is increasing 
monotonically with increasing thickness, y. Based on this isolated result 
lift is maximized for additional thickness of the form assumed by moving 
the position of maximum additional thickness as far aft as allowed and 
introducing as much additional thickness as allowed. 
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Moment Minimization 

„ Minimization o£ the moment coefficient has form 

0 a Min J 

where 

' i ram ji'*i nation resulted in a solution directly opposed to lift 
!!■ sat The position of maximum thickness moved to the forward 
ii> ' iind the amount of additional thickness was minimized. That is, 
the basic NACA 64,-212 airfoil has less adverse moment than any airfoil 

generated by addition of biquadratic thickness to the upper surface of 
the airfoil. 

Lift vs. Moment Trade 

Preliminary work using other airfoil thickness representations has 
indicated that the requirements of lift maximization and moment minimization 
oppose each other. This is confirmed by the results reported above. 

It has been found as a result of previous studies that it is very difficult 

J{ 

to produce an airfoil for which 

C m < .177 - .22 C L 

This function has been used to define airfoil which have favorable lift/ 
moment characteristics by solution of the problem 

f ' 0 = Min [.177 - .22 C L - C M j 

Solution of this problem by directed random-ray and pattern search 
indicates that additional thickness should be added as far forward as 
possible and that maximum amount of additional thickness should be employed. 

Optimisation Summary 

Three optimal airfoil results have been obtained consistent with the 
class of airfoil profiles considered here. These results are summarized 
in Table II. It can be seen that in all cases the position of maximum 
thickness, x, is either at the extreme forward or rearward position allowed. 
Similarly, depending on problem specification, the amount of additional' - 

f 

t 
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rc (i | - <V (x - k) Ap(x)dx 


t ...aicSS is either minimized cn maximized. The low dimensionality of 
this problem (two parameters, x and y) permit a ready mapping of these 
results as a function of x and y. This is done in the following section. 

SYSTEMATIC VARIATION OF 
AIRFOIL SHAPING PARAMETERS 


A systematic investigation on the effect of variations in the airfoil 
shaping parameters x and y was undertaken. The resulting airfoils and 
calculated pressure distributions are presented in Figures 3(a) to 3(v). 

It should bo noted that the airfoils are not drawn to scale in Figure 3. 

To emphasize profile characteristics the vertical scale is exagerated. 

The pressure signatures vary in a radical manner with x and y. The basic 
airfoil exhibits a sharp pressure peak at the leading edge. The magnitude 
of the peak pressure is reduced by introducing additional thickness in 
a forward location, x = .1, and the peak position moves aft. However, if 
the amount of additional thickness is increased the pressure peak magnitude 
again increases. This peak is well aft of the leading edge. This effect 
persists until rearward locations of x are encountered. For example, 
introducing additional thickness at x = .8 results in a rearward "hump" 


in the pressure distribution. The increased circulation produced by this 
hump results in an increased leading edge peak in the airfoil pressure 
distribution. Flow separation would probably be encountered with these 
rearward additional thickness distributions unless devices such as 
rotating cylinders or blowing were employed. 

Figure 4 illustrates the effect of varying position of maximum 
thickness and maximum thickness on lift coefficient. It can be seen that 
lift coefficient is maximized by increasing both x and y. This confirms 
optimization studies in the previous section. Since the additional thick- 
ness and the basic 12% airfoil thickness are additive Figure 3 presents 
lift coefficient as a function of thickness. To first-order the airfoil 
thickness required is 

t/c = 12% + y 
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As it moves to the extremes o£ the range the actual airfoil thickness 
is less than this amount as the positions of maximum thickness on the asic 
and additional thickness distributions are significantly different. 

Moment coefficient variation with x and y is presented in Figure 5. 

It can be soon that the increased lift available from additional thickness 
is accompanied by an increase in undesirable pitching moment coefficient, 
the conclusion of the previous section that moment coefficient is 
minimized by moving x forward and diminishing y is borne out by Figure 5, 
again confirming the optimization study results. 

A final verification of the optimization procedures employed is 
provided by Figure 6. Here the variation of and with x and y is 
presented together with the line function 

.177 - .22 C L - C M = 0 

It can be seen that based on this function the most favorable 

trade involves moving x forward and introducing the maximum y. 

Figure 7 presents the relationship between pressure peak and lift 
coefficient for a range of x and y values. For each value of y (maximum 

V ""i 

additional thickness) there is a point at which the pressure peak magnitude 

is minimized. Cross plotting the peak pressures as a function of C L in 

Figure 8 reveals the minimum peak pressures as a function of C L> 

Figure 9 plots the position of maximum additional thickness as a 

function of C, . As C. increases x moves aft. The associated values of 

y required for the low peak pressure is also plotted in Figure 9. Finally^ 

Figure 10 plots the minimum C attainable as a function of C. using the 

p L. 

biquadratic additional thickness airfoil model. 

n 

CONCLUSION 


A numerical investigation into a class of modified airfoil shapes 
has been completed using full two-dimensional flow potential flow equations. 
Airfoils studied were obtained by modifying- the NACA 64^ -212 airfoil by 
additional thickness distributions based on a biquadratic variation with 
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chordwisc position. Free streamwise Mach number was held constant at 
M « 0.2 and the basic airfoil is held at 6° anglo-of-attack. Results 
of the study may bo summarized as follows: 

1. Significant changos in pressure distribution, lift and 
pitching momont can bo introduced by the biquadratic 
thickness modification. 

2. The requirements for improving lift and moment co- 
efficient characteristics are directly opposed to each 
other. That is, increases in lift result in increases 

in adverse moment. Conversely, decrcnses in adverse moment 
produce decreases in lift. 

3. High lift airfoils require the addition of a thickness 
distribution biased to the rear of the foil and as much 
thickness addition as possible. 

4. Low adverse moments require a thickness distribution 
biased to the front of the foil and as little additional 
thickness as possible. Therefore, the best airfoil based 
on momont considerations is the unmodified foil, 

5. Favorable lift/moment trade-off characteristics are obtained 
by a thickness distribution biased to the front of the foil 
employing as much thickness as possible. 

6. There exists a class of airfoil exhibiting low peak pressures 
for a given which require an intermediate location of 
maximum additional thickness and thickness amount. Generally, 
the position of maximum additional thickness moves to the 

rear with increasing C^, and the amount of additional thickness 
required increases with increasing C L » 
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TABLE I 


CONVERGENCE FOR C L MAXIMIZATION 


M 

C 

jjj 

X 

alfiiac 

y 

n alphai 2) 

r.OOOOE-03 

- C L 
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0 

I 

l 

.2500 r 

-.6941 

1° 

1 

2 
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1‘. 0365003 

-.8942 *' 

10 

1 

3 

,2583 

f. J667O03 

-.6952 

To 

1 

5 

,2626 

1’ 4 155E-03 

-.8968 

2 

1 

6 

,2753 

1,831 IE-03 

-.8982 

2 

1 

7 

* .2879 

2.2466E-03 

-.9022 

2 

t 

3 

.3132 

3 , 0 7 7 6 E « 0 3 

-.9076 

2 

1 

9 

.3633 — 

4',739BE-03 

-.9196-*“ 

1° 

2 

10 

.3688 

4’, 7533E-03 | 

-.9198 

10 

2 

11 

f 2756 

5. 2172E-Q3 I 

-.9216 

10 

2 

12 

. 3824 

5, 561 1 E- 03 ; 

-.9258 

2 

2 

13 

,4010 

6'. 3624 E- 03 

-.9320 

2 

2 

14 
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7 . 2 0 3 7 E - 0 3 

-.9385 

2 

2 

15 

.4567 

8'.6462E-03 * 

-,9520 **“ 


2 

16- 

.5311“ 

1.2131E-02 

-.9826 

10 

3 

17 

,5505 

1.2761E-02 

-.9863 ' 

lo 

3 

16 

.5675 

1 ‘,28 1 2E-02 

-.9920 

10 

3 

20 

,5806 

1.349/IE-02 . 

-.9993 

2 

3 

21 

.6305 

l’,4656E-02 

-1,016 

a 

3 

22 
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• T.6219E-02 

-1,041 

2 

3 

23 
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-1.105 

2 

3 

24 
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* 2‘.439/IE-02 

-1.303 ' 

10 

4 

25 
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2.501 6 E~ 02 

-1,313 
-1,350 - 

i° 

n 

26 
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10 

4 

27 
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3.01 16E-02 

-1.399 

2 

n 
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4 
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Problem 
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64 (I) -212 AIRFOIL AND PRESSURE DISTRIBUTION 



FIGURE 1 


'e) x = .9 


FIGURE 2 
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FIGURE 3(n) 






64(11-062 EXTENSION X-,-6 Y--.06 



FIGURE 3(r) 
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